Magnetic excitations and phonons in the spin-chain compound NaCu202 
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We report an inelastic light scattering study of single-crystalline NaCu2C>2, a spin-chain compound 
known to exhibit a phase with helical magnetic order at low temperatures. Phonon excitations were 
studied as a function of temperature and light polarization, and the phonon frequencies are compared 
to the results of ab-initio lattice dynamical calculations, which are also reported here. The good 
agreement between the observed and calculated modes allows an assignment of the phonon eigen- 
vectors. Two distinct high-energy two-magnon features as well as a sharp low-energy one-magnon 
peak were also observed. These features are discussed in terms of the magnon modes expected in 
a helically ordered state. Their polarization dependence provides evidence of substantial exchange 
interactions between two closely spaced spin chains within a unit cell. At high temperatures, the 
spectral features attributable to magnetic excitations are replaced by a broad, quasielastic mode 
due to overdamped spin excitations. 

PACS numbers: 75.40.Gb, 63.20.-e, 78.30.-j, 75.50.-y 



I. INTRODUCTION 

A large research effort over the past three decades 
has been devoted to the exploration of quantum ef- 
fects in spin-chain materials. At low temperatures, 
most quasi-lD magnets exhibit either gapped spin-liquid 
ground states, or states with collinear long-range mag- 
netic order. 1 In some spin-chain compounds, spin-lattice 
interactions stabilize spin-Peierls ground states. Re- 
cently, several quasi-lD magnets with noncollinear, in- 
commensurate long-range order at low temperatures have 
also been found. LiCu2 02, an insulating compound 
containing edge-sharing Cu02 chains, is a prominent 
example. 2 ' 3,4 ' 5 ' 6 ' 7 Initially, both lattice disorder 3 and spin 
frustration 4 had been discussed as the driving mechanism 
of the observed helimagnetic order. This controversy 
now appears to have been resolved: 6 The incommensu- 
rate magnetic ordering is due to long-range competing 
exchange interactions along the Cu02 chains. Questions 
remain, however, about the microscopic role of the sub- 
stantial number of Li impurities substituting copper ions 
on the spin chains, and about the magnitude of the ex- 
change interactions between two closely spaced Cu02 
chains within a unit cell. Inelastic neutron scattering 
measurements 5 suggest that these interactions are sub- 
stantial, whereas electronic structure calculations 4 indi- 
cate much smaller values. 

Investigations of the disorder-free, isostructural com- 
pound NaCu202 8 have now begun to help resolve these 
and other questions. 9 This material is also interesting be- 
cause it is the end member of the series Nai +2; Cu202. For 
nonzero x, the spin chains of these materials are doped 



and exhibit intriguing phenomena such as Wigner crys- 
tallization of holes. 10 Like its Li-analog, NaCu202 has an 
orthorhombic crystal structure with space group Pmna. 
The lattice parameters are a= 6.2087 A, b=2.9343 A, and 
c=13.0548 A. 8 Magnetic Cu 2+ ions form two inequivalent 
edge-sharing copper oxide chains per unit cell. They run 
along the b axis and are shifted relative to each other 
by 6/2. Along the a and c-directions, they are sepa- 
rated by nonmagnetic Na + and Cu + ions, respectively. 
The Cu-O-Cu bond angle is 92.9°. This results in a 
small, ferromagnetic, nearest-neighbor (NN) exchange in- 
teraction. The antiferromagnetic next nearest-neighbor 
(NNN) interaction is the largest parameter in the spin 
Hamiltonian. An analysis of high-temperature suscepti- 
bility demonstrates the significance of interactions reach- 
ing up to fourth-nearest neighbors of a Cu 2+ ion. 9 The 
static susceptibility displays a broad maximum around 
52 K, associated with the evolution of short-range anti- 
ferromagnetic correlations. 9 Upon cooling, susceptibility 
and specific heat measurements show two magnetic phase 
transitions at 12 and 8 K. Neutron diffraction results 9 
demonstrate that the former transition is due to a non- 
collinear magnetic ordering with the propagation vector 
(0.5, C, 0) with C = 0.227, similar to that of LiCu 2 2 
where £ = 0.174. 3 This underscores the conclusion that 
lattice disorder (which is negligible in NaCu202) is not 
responsible for the helimagnetism in this class of com- 
pounds. The ordered moment of 0.56 fiB is reduced with 
respect to the full spin-1/2 moment due to quantum fluc- 
tuations. The latter transition might be related to spin 
canting. 

While both Raman scattering 7 and inelastic neutron 
scattering 5 work on magnetic excitations in LiCu202 has 
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been reported, the magnetic dynamics of NaCu2C>2 has 
thus far not been investigated. In this article, we report 
an investigation of magnetic excitations in NaCu2 02 by 
polarized Raman scattering. In the magnetically long- 
range ordered state, three distinct features attributable 
to magnetic excitations are observed. Two modes around 
18 and 58 cm -1 are observed both with parallel in-chain 
and crossed polarizations. This points to significant in- 
terchain interactions. In addition, a broad maximum 
around 390 cm -1 is present in in-chain polarization only. 
The temperature dependences of the low-energy and the 
two high-energy modes are characteristic of one- and 
two-magnon excitations, respectively. Upon heating, the 
three magnetic modes disappear, and we observe a pro- 
nounced quasielastic response predominantly in the in- 
chain polarization channel. These results are discussed 
in the light of previous work on CuC>2 spin chain materi- 
als, especially recent neutron scattering data and model 
calculations on LiCu202- 5 

The phonon spectrum of NaCu2 02 was also investi- 
gated as a function of temperature and light polariza- 
tion. The observed modes agree well with those pre- 
dicted by ab-initio lattice dynamical calculations, which 
allows us to identify the phonon eigenvectors at the zone 
center. Unusually strong two-phonon modes with an un- 
conventional temperature dependence are also observed. 
These features may be related to a broad electronic mode 
around 500 cm -1 , which possibly arises from electroni- 
cally active defects. 

This article is organized as follows. In Section II, we 
describe technical details pertaining to the Raman exper- 
iments and the lattice dynamical calculations. Sections 
III. A and III.B contain Raman data on phonon and mag- 
netic excitations, respectively, along with a discussion of 
possible mode assignments. Section IV provides a sum- 
mary of our findings. 



II. EXPERIMENTAL DETAILS AND LATTICE 
DYNAMICAL CALCULATIONS 

The single crystals were grown using a self- flux method 
as described elsewhere. 8 The chemical composition of the 
studied sample is close to the ideal stoichiometry, as dis- 
cussed previously. 9 Raman scattering measurements were 
performed in quasi-backscattering geometry with the ex- 
citation line A = 514.5 nm of an Ar + laser. The laser 
power of 7 mW was focused to a 0.1 mm diameter spot 
on the sample surface. Stokes and anti-Stokes measure- 
ments confirmed that the incident radiation did not in- 
crease the temperature of the samples by more than a few 
K. Spectra of the scattered radiation were collected by a 
DILOR-XY triple spectrometer and recorded by a nitro- 
gen cooled charge-coupled device detector with a spectral 
resolution less than ~ 1 cm -1 . 

The optical phonon frequencies were calculated within 
the density functional theory using the VASP software 
package. 11 Projector augmented wave (PAW) potentials 



with an energy cutoff of 283, 259, and 273 eV were used 
for O, Na, and Cu, respectively. The atomic coordinates 
were first optimized in the single unit cell. Then a (2,4,1) 
supercell, which is 2, 4, and 1 times larger in the respec- 
tive a, b and c directions, was constructed and the co- 
ordinates of the 160 atoms were optimized again. For 
the dispersion relation calculations, two k-points in the a 
and b directions and one in the c direction were retained, 
while for the (2,4,1) supercell a single k-point was used. 
Supercell calculations are required in order for the mag- 
nitude of the force constants to become negligible (which 
happens beyond 10 Angstroms) as well as to calculate 
the phonon dispersions exactly instead of using an in- 
terpolation. A series of single point energy calculations, 
which give the Hellmann-Feynmann forces acting on all 
atoms in the supercell, were then performed on the struc- 
ture obtained by displacing, one at a time, each of the 
20 inequivalent atoms in positive and negative directions 
along the Cartesian axes. Finally, the dynamical matrix 
for any point in reciprocal space was generated and di- 
agonaliscd. 



III. RESULTS AND DISCUSSION 
A. Phonons 

Figure 1 displays Raman spectra in bb-, ba-, and aa- 
polarizations at 5 and 295 K. Most of the modes shown 
arc of lattice vibrational origin. In order to assign these 
modes, we performed a group theoretical analysis as well 
as ab-initio lattice dynamical calculations following the 
description in Section II. NaCu202 has the space group 
Pnma, where all atoms have a 4c site symmetry. Apart 
from the acoustic (Bi„ + B2„ + B3„) modes, the factor 
group analysis yields the following Raman- and infrared- 
active modes: 10 A 9 (aa,bb,cc) + 5 B lg (ab) + 10 B 2s (ac) 
+ 5 B 3s (bc) + 5 A„ + 9 Bi„ + 4 B 2u + 9 B 3u . Thus, 
in bb- and aa-geometries we expect 10 A g modes, and in 
the ba- configuration 5 Bi g modes. 

The observed peak frequencies are listed in Table I 
together with the calculated ones. Given that the calcu- 
lation was performed without adjustable parameters, the 
agreement of observed and calculated frequencies is quite 
satisfactory. We can thus tentatively assign the observed 
modes to specific vibration patterns. Figures 2 and 3 pro- 
vide a synopsis of the eigenvectors of the Raman-active 
modes resulting from the calculation. It can be seen that 
the highest-energy modes around 500 cm -1 are predomi- 
nantly oxygen vibrations, as expected based on the small 
mass of the oxygen ion. 

In both aa- and ^^-polarizations we observe a total 
of 13 A g modes. This number is larger than the ex- 
pected 10 modes. Specifically, the highest-energy oxygen 
vibration around 560 cm -1 seems to be split into three 
modes, indicating different environments of oxygen ions 
in the crystal structure (see Fig. 7). This is either due 
to coherent deviations from the reported Pnma lattice 
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TABLE I: Comparison between experimental and calculated 
phonon frequencies of the Raman-active modes in a unit of 
cm -1 . 
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symmetry too subtle to be discernible in the reported 
x-ray and neutron 9 powder diffraction data, or to lat- 
tice defects. In afe-polarization 11 modes are observed, 
a number that is also larger than the symmetry-allowed 
5 Bi s modes. Our calculations show that the estimated 
frequency separation between the A g and B iff modes is 
only 1 — 2cm~ 1 . This lies within our spectral resolution 
of <~ 1 cm -1 . It is therefore difficult to assess whether the 
modes observed in a6-polarization are pure Bi 3 modes, or 
whether they are due to an admixture of A g modes due 
to polarization leakage caused by lattice defects and/or 
other mechanism. 

The temperature dependence of the peak frequencies 
and linewidths of representative phonons, extracted by 
fitting the corresponding spectra to Lorentzian profiles, is 
shown in Fig. 4. Nearly all phonon modes exhibit a sim- 
ilar behavior. With decreasing temperature the peak fre- 
quencies harden and the linewidths narrow gradually, be- 
fore saturating at about 70 K. This behavior is commonly 
observed as a consequence of lattice anharmonicity. How- 
ever, two phonons exhibit phenomena not attributable to 
anharmonic effects alone. First, the lincwidth of the 146 
cm -1 mode, which according to the calculation mostly 
involves a 6-axis vibration of the Cu 1+ ions against the 
rest of the lattice (Fig. 3), exhibits a subtle but distinct 
drop upon cooling from 50 K to the lowest temperatures 
(see the arrow in Fig. 4). As the Cu02 chains vibrate in 
unison, a strong spin-phonon coupling is not expected, 
and the relationship of this subtle anomaly to the devel- 
opment of magnetic correlations is not obvious. Second, 
the 502 cm -1 oxygen vibration shows an asymmetric line- 
shape, which appears to be due to an anti-resonance with 
a broad electronic mode centered around 500 cm -1 . The 
origin of the electronic mode will be addressed in Section 
III.B below. 



Here we draw attention to the unusual behavior of the 
two-phonon peak around 1000 cm -1 , which is probably 
an overtone of the 502 cm -1 mode (Fig. 1). In 66- 
polarization and at 5 K, the integrated intensity of the 
two-phonon peak is nearly an order of magnitude higher 
than that of the one-phonon peak. At room temperature, 
the two-phonon peak decreases strongly in intensity. In 
contrast, the ratio of the integrated intensities of one- 
and two-phonon modes for the aa-configuration is only 
~ 0.3 in the temperature regime investigated. Since two- 
phonon excitations are usually weak and hardly change 
with increasing temperature, the behavior of this mode 
is quite surprising. Qualitatively similar (albeit less pro- 
nounced) effects were observed in corner-sharing cop- 
per oxide chain systems, and interpreted in terms of 
a Frohlich interaction mechanism. 12,13 This mechanism 
can also lift the selection rules such that infrared-active 
modes may become observable by Raman spectroscopy. 
This might explain the observed extra phonon modes. A 
quantitative understanding of the two-phonon modes of 
NaCu 2 2 remains a subject of future work. 



B. Magnetic excitations 

Magnetic excitations in the Raman spectra can be 
identified by virtue of their temperature and polarization 
dependences. The upper panel of Fig. 1 shows a sub- 
stantial temperature-induced rearrangement of spectral 
weight in in-chain 66-polarization. The low-energy part 
of the spectrum that is affected by this redistribution is 
highlighted in Fig. 5. Upon cooling, a broad quasielastic 
background is progressively replaced by a more struc- 
tured spectrum. At 5 K, three features attributable to 
magnetic excitations emerge: (i) an ~ 50 cm -1 wide peak 
centered at 390 cm -1 ; (ii) a weaker peak with an asym- 
metric lineshape extending from about 50 to 170 cm -1 ; 
(iii) and a sharp peak around 17 cm -1 . The latter two 
features are also present in a6-polarization (see the mid- 
dle panel of Fig. 1). Their temperature dependence is 
highlighted in Figs. 6 and 7, respectively. The magnetic 
Raman operator is given by 1Z oc Si-Sj, 1 so that magnetic 
Raman scattering is dominant in the polarization where 
there are substantial exchange paths. The observation of 
features (ii) and (iii) in both 66- and 6a-polarizations thus 
indicates that they are appreciably affected by interchain 
coupling. 

Before discussing these features in detail, we recall the 
salient properties of the spin Hamiltonian of NaCu 2 2 . 9 
The largest exchange coupling parameter is the antiferro- 
magnetic NNN coupling along the spin chains ( J 2 = 62.6 
cm -1 ). The helix structure is induced by a frustrating 
ferromagnetic NN interaction J\ = 11.4 cm -1 . Longer- 
range exchange interactions along the chains and/or in- 
teractions between spins within a closely spaced pair of 
spin chains are also required to obtain a good fit to the 
magnetic susceptibility and the propagation vector of the 
helix structure. 
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The magnetic excitations of NaCu202 have thus far 
not been studied, but an inelastic neutron scattering 
study has recently addressed the low-energy spin excita- 
tions in LICU2O2. 5 Based on model calculations, 14 the 
largest parameter in the spin Hamiltonian, J2, is ex- 
pected to be very similar in the two systems. From 
the ~30% smaller pitch angle of the magnetic helix in 
NaCu2C>2 one can infer a correspondingly smaller NN 
coupling J\. This is expected, because the Cu-O-Cu 
bond angle of NaCu 2 2 is closer to 90°, and it is con- 
sistent with the data currently available on this com- 
pound. Although there may be further subtle differences 
between the spin Hamiltonians of both compounds, the 
neutron scattering data and model calculations reported 
for LiCu2C>2 can thus serve as a useful guideline for the 
interpretation of our data. 5 Note that the Raman data 
presented here cover the high-energy and low-energy seg- 
ments of the spin wave spectrum not probed by the neu- 
tron experiments on LiCu 2 2 . Both data sets can thus 
be regarded as complementary in this respect as well. 

We first discuss the two broad modes centered at 70 
and 390 cm -1 . The temperature dependence of the for- 
mer feature is shown in more detail in Fig. 5. It per- 
sists well into the paramagnetic state and merges into 
the background at temperatures exceeding about 80 K. 
The 390 cm -1 feature persists to even higher tempera- 
tures in excess of 150 K, without a strong renormaliza- 
tion of its lineshape (Fig. 7). This is the temperature 
dependence expected for two-magnon Raman scattering 
in low-dimensional magnets, where short-range spin cor- 
relations persist well into the paramagnetic state. The 
presence of the higher-energy two-magnon peak in 66- 
polarization, and its absence in a6-polarization, confirm 
that the magnetic zone-boundary energy is determined 
predominantly by exchange interactions within a single 
spin chain, as expected based on the discussion above. 
Its energy is expected to be slightly below twice the mag- 
netic zone-boundary energy, due to magnon-magnon in- 
teractions. This is indeed the the zone-boundary 
energy computed in Rcf. 5 is 26 meV (=210 cm -1 ). 

The presence of a second two-magnon feature is more 
unusual, but it can be understood as a consequence of 
the helical spin correlations present in NaCu202 and 
LiCu2C>2. The helical spin structure and the four- atom 
basis in the magnetic unit cell give rise to a complex 
spin wave spectrum (Fig. 6a of Ref. 5), with several 
magnon branches exhibiting cither maxima or minima in 
the energy range 8-12 meV (65-97 cm -1 ) for LiC^C^- 
Because of the lower value of J\ (see discussion above), 
this energy scale is expected to be reduced in NaCu 2 02- 
The high two-magnon density of states due to the van- 
ishing spin-wave velocity at these extrema of the disper- 
sion relation naturally explains the second two-magnon 
Raman peak. The observation of this peak in both 66- 
and a6-polarization indicates that magnons in the energy 
range around 10 meV are significantly influenced by in- 
terchain exchange interactions within a pair of closely 
spaced chains. This agrees with the findings of Ref. 



5, but seems incompatible with the prediction of neg- 
ligible interchain exchange interactions based on LDA 
calculations. 4 

We will now address the origin of the peak at 18 cm -1 . 
As shown in detail in Fig. 6, the peak broadens and 
shifts to lower frequency upon heating, and it vanishes 
above the magnetic transition temperature. This behav- 
ior is characteristic of one-magnon excitations. A one- 
magnon excitation with almost identical energy has been 
observed in CuGeC>3 , a spin-Peierls compound also based 
on edge-sharing copper oxide chains, in magnetic fields 
high enough to stabilize a phase with an incommensu- 
rate modulation of both the spin density and the crystal 
lattice. 15 However, the presence of both spin and lattice 
modulations has made it difficult to arrive at a conclusive 
interpretation of the nature of this mode and the origin of 
its Raman-activity. 15 ' 16,17 One-magnon modes of compa- 
rable energies have also been observed in Raman scatter- 
ing experiments on the two-dimensional antiferromagnet 
La2Cu04. 18 These zone-center transverse magnon modes 
are Raman active by virtue of spin-orbit coupling, which 
manifests itself in an exchange anisotropy two orders of 
magnitude smaller than the isotropic Cu-O-Cu superex- 
change interaction. 

The most straightforward interpretation of the 18 
cm -1 peak is based on a transverse magnon at q = 0. 
This is consistent with neutron scattering measurements 
on LiCu202, which tentatively indicate a magnon gap of 
comparable magnitude. 5 If this scenario is valid, the size 
of this gap constitutes a large fraction of the isotropic 
superexchange interaction. This may at first seem sur- 
prising in view of the nearly quenched orbital moment of 
the Cu 2+ ion. However, prior work on edge-sharing cop- 
per oxide chain systems has experimentally demonstrated 
substantial superexchange anisotropics. 19 The surpris- 
ingly large relative magnitude of the non-isotropic part 
of the exchange interaction can be understood based 
on electronic structure calculations for the 90° Cu-O- 
Cu bond geometry. 19 A zone-center transverse magnon 
is therefore a plausible origin of the 18 cm -1 peak. 

An alternative interpretation of this peak is in terms 
of a longitudinal magnon. Such excitations are usu- 
ally heavily damped, but can live long enough to be 
probed near quantum critical points. Indeed, longitudi- 
nal magnons have recently been observed in Raman scat- 
tering experiments 20 on the coupled spin-tetrahedra sys- 
tem Cu2Te205Br2, which exhibits incommensurate heli- 
cal order and a spin excitation spectrum not unlike those 
of NaCu202- 21 ' 22 Because of the sizable ordered moment 
of NaCu202, 9 this scenario is less plausible. However, 
a conclusive assignment of the 18 cm" 1 peak will have 
to await neutron spectroscopy experiments on NaCu202 
akin to those that have led to the identification of longi- 
tudinal magnons in other systems. 23 ' 24 ' 25 Another open 
question concerns the temperature dependence of the fre- 
quency and intensity of the 18 cm -1 mode. In contrast 
to one-magnon excitations in most collinear magnets, 
whose frequency and intensity roughly scale with the or- 
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der parameter, we observe a temperature-linear behavior 
of these quantities (Fig. 6). Detailed model calculations 
for the observed helical magnetic structure are required 
to address the origin of this behavior. 

At high temperatures, the one- and two-magnon peaks 
are replaced by an intense quasiclastic peak due to 
overdamped spin excitations. Similar phenomena have 
been studied in other quasi-one-dimensional magnets 
such as CuGeC>3. 26,27 Using a hydrodynamic descrip- 
tion of the correlation function, 1 the quasiclastic Raman 
peak is well described by a Lorentzian profile I{uo) oc 
C m T 2 D T k 2 /[ui 2 + (D T k 2 ) 2 }, where k is the scattering 
wave vector, Dt the thermal diffusion constant, and C m 
the magnetic specific heat. Figure 8 displays a represen- 
tative fit of the quasielastic response to a Lorentzian pro- 
file at 295 K. The inset of Fig. 8 shows the temperature 
dependence of the scattering amplitude, which increases 
rapidly with increasing temperature and then saturates 
around 200 K. This is the same temperature at which the 
higher-energy two-magnon peak vanishes, heralding the 
disappearance of short-range spin-spin correlations. Here 
we note that deviations from the Lorentzian profile are 
discernible above 100 cm -1 (not shown here). While this 
is partly due to an ambiguity in determining the scatter- 
ing background, it may also reflect interference from an 
additional scattering channel at higher energies, which 
we will now discuss. 

Figure 7 indeed shows a broad, weakly temperature 
independent peak centered around 500 cm -1 , the tail of 
which may well extend down to about 100 cm -1 . The 
peak is observed almost exclusively in the 66-polarization 
channel, although some traces may also be present in 
the other polarizations. As already mentioned in Section 
III. A, the oxygen vibration at 502 cm -1 exhibits a pro- 
nounced anti-resonance behavior with this mode. The 
peak does not shift as the laser frequency is varied, so 
that luminescence can be ruled out. Since the strong 
anisotropy, the nearly temperature dependent spectral 
weight, and the featureless shape of the peak are in- 
consistent with multiphonon scattering, it is probably of 
electronic origin. 

An unambiguous assignment of this peak is not pos- 
sible based on the data at hand. A priori, one has to 
consider orbital, charge, and magnetic excitations of the 
electrons. Raman scattering experiments have recently 
uncovered evidence of orbital excitations in transition 
metal oxides with t2 g valence electrons at energies rang- 
ing from ~ 500 to ~ 2000 cm -1 , 28,29 but orbital excita- 
tions for the e g electron of a Cu 2+ ion are expected to 
occur at much higher energies. Because of the presence 
of both Cu 1+ and Cu 2+ ions in NaCu 2 02, one might also 
consider low-energy charge excitations as the origin of the 
peak. The eigenvectors of both of the phonons exhibit- 



ing anomalous behavior suggest that the corresponding 
vibrations modulate Cu 1+ -Cu 2+ charge-transfer excita- 
tions. However, band structure calculations and photoc- 
mission data for LiCu202 indicate considerably larger en- 
ergies for such excitations. 30 Finally, a spinon continuum 
extending beyond the magnon bandwidth is expected on 
general grounds in one-dimensional magnets. Indeed, 
Raman scattering from such continuum excitations has 
been reported for edge-sharing copper oxide chains. 12 In 
this scenario, the phonon anomalies discussed in Section 
III. A may be consequences of a spin-phonon interaction, 
also invoked to explain phonon anomalies in CuGeC>3. 31 
While we consider this as the most likely interpretation 
of the 500 cm -1 peak, we cannot rule out entirely that 
it is a consequence of defects in the NaCu202 crystal 
structure. As the Na/Cu ratio determined in prior work 
on NaCu202 was stoichiometric within an error mar- 
gin of about 2%, 9 defects on these atomic sites can be 
present only in very small concentrations. Oxygen de- 
fects can, however, not be ruled out. Whatever their 
origin, it is conceivable that electronically active lattice 
defects could give rise to charge carriers in shallow traps, 
which shake off spin excitations propagating along the b- 
direction when excited into the conduction band, giving 
rise to a broad Raman peak. Further work with other 
probes is required to conclusively establish the origin of 
the 500 cm -1 peak. 



IV. CONCLUSIONS 

To conclude, we have presented a detailed study of the 
phonon spectrum and the magnetic response of the heli- 
cally ordered spin-chain system NaCu202 using Raman 
spectroscopy. Based on lattice dynamical calculations, 
most of the phonon modes could be assigned. Several 
spectral features originating from magnetic excitations 
were observed as well. The rich magnetic Raman spec- 
trum of NaCu 2 2 is unusual for a quasi-one-dimensional 
magnet. However, a comparison to model calculations 
and neutron scattering data on the isostructural com- 
pound LiCu2C>2 and prior Raman scattering work on 
other copper oxide chain systems yields plausible inter- 
pretations of all of the observed modes. 
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FIG. 1: (Color online) Polarization dependence of Raman 
spectra at 5 and 295 K in bb-, ba-, and aa- polarizations, 
respectively. The downward vertical arrows indicate magnetic 
excitations while the upward one denotes a Fano-shaped dip 
corresponding to the phonon at 502 cm -1 . 
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FIG. 2: (Color online) Calculated normal frequencies and 
eigenvectors of the ylg-symmetry modes. The relative am- 
plitude of the vibrations (arrow length) is correct within each 
mode. For the sake of clarity a multiplication factor has been 
applied to some modes. The red balls stand for the oxygen 
ions, the blue ones for sodium, the yellow ones for the Cu +2 , 
and the orange ones for the Cu +1 . 




FIG. 3: (Color online) Calculated normal frequencies and 
eigenvectors of the Bi 3 -symmetry modes. The relative am- 
plitude of the vibrations (arrow length) is correct within each 
mode. For the sake of clarity a multiplication factor has been 
applied to some modes. 
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FIG. 4: (Color online) (Left panel) Temperature dependence 
of the peak frequencies of selected phonons. (Right panel) 
Phonon linewidths as a function of temperature. 
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FIG. 5: (Color online) Temperature dependence of magnetic 
Raman spectra in ba-polarization at 5, 20, 40, 60, and 80 K 
(from upper to lower curve). 
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FIG. 6: (Color online) (a) Low-frequency magnetic Raman 
spectra in 6a-polarization as a function of temperature. (b,c) 
Temperature dependence of the peak position and integrated 
intensity on a reduced temperature scale, respectively. 




FIG. 7: (Color online) Raman spectra in 66-polarization as a 
function of temperature; 290, 150, 80, 50, and 5 K from upper 
to lower curve. 
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FIG. 8: (Color online) Result of a fit of the quasielastic re- 
sponse to a Lorentzian profile at 295 K. Inset: Temperature 
dependence of the integrated scattering intensity of this pro- 
file. 



